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Abstract: We present a thorough analysis for the quasinormal (QN) behavior, as- 
sociated with the decay of scalar, electromagnetic and gravitational perturbations, of 
Schwarzschild black holes in anti-de Sitter (AdS) spacetimes. As it is known the AdS 
QN spectrum crucially depends on the relative size of the black hole to the AdS radius. 
There are three different types of behavior depending on whether the black hole is large, 
intermediate, or small. The results of previous works, concerning lower overtones for 
large black holes, are completed here by obtaining higher overtones for all the three black 
hole regimes. There are two major conclusions that one can draw from this work: First, 
asymptotically for high overtones, all the modes are evenly spaced, and this holds for 
all three types of regime, large, intermediate and small black holes, independently of 
I, where I is the quantum number characterizing the angular distribution; Second, the 
spacing between modes is apparently universal, in that it does not depend on the field, 
i.e., scalar, electromagnetic and gravitational QN modes all have the same spacing for 
high overtones. We are also able to prove why scalar and gravitational perturbations are 
isospectral, asymptotically for high overtones, by introducing appropriate superpartner 
potentials. 
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1 Introduction 



Any physical system has its modes of vibrations. For non-dissipative systems these modes, 
forming a complete set, are called normal. Each mode having a given real frequency of 
oscillation and being independent of any other. The system once disturbed continues to 
vibrate in one or several of the normal modes. On the other hand, a black hole, as any 
other gravitational system, is a dissipative system since it emits gravitational radiation. 
One has to consider, instead, quasinormal modes (QNMs) for which the frequencies are 
no longer purely real, showing that the system is loosing energy. QNMs are in general not 
complete and though insufficient to fully describe the dynamics, contain a great amount 
of information. For instance, they dominate the signal during the intermediate stages 
of the perturbation. Indeed, calculations ranging from the formation of a black hole in 
gravitational collapse to the collision of two black holes provide clear evidence that no 
matter how one perturbs a black hole, its response will be dominated by the QNMs. 
Through the QNMs, one can also probe the black hole mass, electric charge and angular 
momentum by inspection of their characteristic waveform, as well as test the stability of 
the event horizon against small perturbations. Moreover, the interest in QNMs has now 
broaden, they might be related to fundamental physics. In the context of black holes in 
asymptotically flat spacetimes the importance of QNMs has long ago been recognized. 
Due to a crescent increase of interest for black holes in asymptotically de Sitter (dS) and 
asymptotically anti-de Sitter (AdS) spacetimes the study of QNMs has now spread into 
these. 

In asymptotically flat spacetimes the idea of QNMs started with the work of Regge 
and Wheeler [Q where the stability of a black hole was tested, and were actually first nu- 
merically computed by Chandrasekhar and Detweiler several years later |2j. It continues 
to be a very active field due to the eminent possibility of detecting gravitational waves 
from astrophysical sources. There are two standard reviews in the field jSj- QNMs in 
asymptotically flat spacetimes have recently acquired a further importance since it has 
been proposed that the Barbero-Immirzi parameter, a factor introduced by hand in order 
that Loop Quantum Gravity reproduces correctly the black hole entropy, is equal to the 
real part of the quasinormal (QN) frequencies with a large imaginary part jlj (see [5] for 
a short review). For further developments in calculating QN frequencies in Kerr space- 
times and in asymptotically flat black holes spacetimes in d-dimensions see [Bj. Some 
other recent calculations of QN frequencies in asymptotically flat black hole spacetimes 
can be found in jjj. 

In asymptotically dS spacetimes the calculation of QN frequencies was first done by 
Moss and collaborators [H] in which the stability of the Cauchy horizon of a charge black 
hole was analyzed. In Cardoso and Lemos jH] an analytical method was devised to study 
the case in which the black hole and the cosmological horizons are very close to each other. 
This analytical method has recently been extended to higher dimensional Schwarzschild 
dS black holes 10 and also to higher order in the difference between the cosmological 
and horizon radius [TT] . A different analytical approach has also been developed in [T2j . 

In asymptotically AdS spacetimes, which is the background spacetime for our paper, 
there has been a great amount of work because the AdS/CFT correspondence conjecture 
[T3] makes the investigation of QNMs important. According to it, the black hole cor- 
responds to a thermal state in the conformal field theory (CFT), and the decay of the 
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test field in the black hole spacetime, corresponds to the decay of the perturbed state 
in the CFT. The dynamical timescale for the return to thermal equilibrium can be done 
in AdS spacetime, and then translated onto the CFT, using the AdS/CFT correspon- 
dence. Many authors have now delved into these calculations in several different black 
hole settings in several different dimensions (see [14] for a sample). In this paper we are 
interested in the 4-dimensional Schwarzschild-AdS spacetime. The lowest lying modes 
(i.e., the less damped ones parametrized by the overtone number n — 0) for this space- 
time were found by Horowitz and Hubeny and completed by Konoplya JH] for the 
scalar case, and by Cardoso and Lemos for the electromagnetic and gravitational 
case. Recently, Berti and Kokkotas have confirmed all these results and extended 
them to Reissner-Nordstrom-AdS black holes. Here, we shall take a step further in car- 
rying on this program by computing numerically, through an extensive search, the high 
overtone QN frequencies for scalar, electromagnetic and gravitational perturbations in 
the Schwarzschild-AdS black hole. We shall do an extensive search for the high overtone 
QN frequencies, (n > 1). We find that the modes are evenly spaced for frequencies with 
a large imaginary part. Moreover the scalar, electromagnetic and gravitational pertur- 
bations all possess, asymptotically for high overtones n, QN frequencies with the same 
spacing, and this spacing is /-independent. While we can numerically prove this with 
great accuracy for large black holes, it remains just a conjecture for small and intermedi- 
ate black holes. We shall also see that the QN frequencies of the toroidal black hole with 
non-trivial topology ^Hj are identical to the QN frequencies of a large Schwarzschild-AdS 
black hole (20]. 



2 Equations and Numerical Method 

We shall deal with the free evolution of massless classical fields in the background of a 
Schwarzschild-AdS spacetime with metric given by 

dv 2 

ds 2 = f(r)dt 2 - ^— - r 2 (d9 2 + sin 2 6d<f) 2 ) , (1) 
f( r ) 

where, f(r) = jg + l — ™,Ris the AdS radius, and M the black hole mass (Newton's 
constant Gn and the velocity of light are set to one). The evolution of scalar, electromag- 
netic, and gravitational fields can be followed through the Klein-Gordon, Maxwell and 
Einstein equations, respectively. If one assumes that the fields are a small perturbation 
in the background given by equation (jTJ). then all the covariant derivatives can be taken 
with respect to the metric (JTJ) . It is then possible to show that the evolution equations 
are all of the same type, i.e., a second order radial differential equation (for more details 
we refer the reader to |T3] for the scalar case and to J7j for the electromagnetic and 
gravitational case). The wave equation is 



+ 



drl 

where the tortoise coordinate r* is defined as 

dr 



uo 2 - V(r) *(r) = 0, (2) 
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and the potential V appearing in (j! 
Explicitly, for scalar perturbations, 

V s = f(r) 



depends on the specific field under consideration. 



1(1 + 1) 



2M 



+ 



R 2 



while for electromagnetic perturbations, 



fir) 



1(1 + 1) 



(4) 



(5) 



The gravitational perturbations decompose into two sets [T7j, the odd and the even parity 
one. For odd perturbations the potential V(r) in (J2J) is 



V 



odd 



/(*") 



1(1 + 1) 6M 



(6) 



while for even perturbations, we have 

2 /( r ) 9M 3 + 3a 2 Mr 2 + a 2 (1 + a) r 3 + 3M 2 (3ar + 3-J 



(3M + ar) 



(7) 



where a — ~ [/(/ + 1) — 2]. In all cases, we denote by I the angular quantum number, that 
gives the multipolarity of the field. We can of course rescale r, r — > ^. If we do this, the 
wave equation takes again the form (J2J) with rescaled constants i.e., r + — > to — > uR, 
where r + is the horizon radius. So, we can take R = 1 and measure everything in terms 
of R, the AdS radius. Eq. (J2J should be solved under appropriate boundary conditions, 
i.e., incoming waves near the horizon, 



+ > 



and no waves at infinity, 



* = 



-> oo . (9) 

We note that there are other reasonable boundary conditions at infinity, in particular for 
the gravitational perturbations. For instance, one can define Robin boundary conditions 
in such a way as to preserve certain dualities between the odd and the even gravita- 
tional perturbations. However, it was verified numerically by Moss and Norman that 
Dirichlet or Robin boundary conditions yield approximately the same result, so we shall 
keep \I/ = , r — >• oo. Moreover, Cardoso and Lemos J7j proved that for high overtone 
QN frequencies the duality is preserved, so in this regime the distinction is irrelevant. 
Thus, to compute the QN frequencies to such that the boundary conditions (jHJ) and © 
are preserved, we follow the Horowitz- Hubeny approach ^j. Within this approach we 
need to expand the solution to the wave equation around x + = ^- (x = 1/r), 



= J2 a k(u)(x-x+y 

k=0 



(10) 



and to find the roots of the equation ty(x = 0) = 0. First, one should substitute fTUj) 
into the wave equation in order to obtain a recursion relation for ^H]- Then, one 



5 



has to truncate the sum (|10p at some large k = N and check that for greater k the roots 
converge to some true root which is the sought QN frequency. The higher the overtone 
number, and the smaller the black hole size, the larger the number N at which the roots 
of the equation = 0) = converge. Yet, since in the series (jlUj) each next term 
depends on all the preceding terms through the recursion relations, when N is too large, 
the tiny numerical errors in the first terms, start growing as N ~ 10 2 — 10 3 or greater. As 
a result the roots suffer a sharp change for a small change on any of the input parameters, 
displaying a "noisy" dependence. To avoid this we have to increase the precision of all 
the input data and the recursion relation we are dealing with from the standard 20-digit 
precision up to a precision such that further increasing of it will not influence the result 
for the QN frequency. For small black holes the roots start converging at very large N 
only, for instance, when r + = 1/20 we can truncate the series at N ~ 3 • 10 4 , but not 
before. Since for finding roots of © we have to resort to trial and error method, the 
above procedure consumes much time, but nevertheless allows to compute QNMs of small 
AdS black holes and to find the higher overtones we are seeking. 

3 Numerical Results 

In this section we will present the numerical results obtained using the numerical pro- 
cedure just outlined in the previous section. The results will be organized into three 
subsections: scalar, electromagnetic and gravitational perturbations. For each field, we 
shall also divide the results into three different regimes: large, intermediate and small 
black holes, since the results depend crucially on the regime one is dealing with. Here a 
large black hole stands for a black hole with r + ^> 1, an intermediate black hole is one 
with r + ~ 1, and a small black hole has a horizon radius r + < 1. We shall then try to 
unify these results. For each horizon radius r + and angular quantum number I there is an 
infinity of QN frequencies (or overtones). We shall order them according to the standard 
procedure, by increasing imaginary part. Accordingly, the fundamental QN frequency 
is defined as the one having the lowest imaginary part (in absolute value) and will be 
labeled with the integer n = 0. The first overtone has the second lowest imaginary part 
and is labeled with n = 1, and so on. The QN frequencies also have a real part, which 
in general display an increase along with the imaginary part. To the lowest value of the 
imaginary part corresponds the lowest value of the real part, to the second lowest value 
of the imaginary part corresponds the second lowest value of the real part, and so on. 
Thus n, the overtone number, is also a number that in general increases with the real 
part of the frequency (or energy) of the mode. This seems to be a characteristic of AdS 
space only, due to the special boundary conditions associated with this spacetime. This, 
in a sense, is to be expected since the wave equation to be studied is a Schrodinger type 
equation, where for quantum non-dissipative bound systems, such as the hydrogen atom 
or a particle in an infinite well potential, the principal quantum number n (which here 
has been called the overtone number) appears due to the boundary conditions of the ra- 
dial equation, a typical eigenvalue problem, and is related directly with the frequency of 
vibration of the orbital. The similarity is not full, though, since the boundary condition 
at the black hole is of a different kind. However, for pure AdS spacetimes, when there is 
no black hole and the boundary conditions are of infinite well type, the overtone number 
n is is indeed a principal quantum number (see Appendix A). 
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3.1 Scalar quasinormal frequencies 

The fundamental scalar QN frequencies were first computed by Horowitz and Hubeny 
[To] for intermediate and large black holes. Konoplya JH] extended these calculations 
to the case of small black holes. Recently Berti and Kokkotas ^H] rederived all these 
results. Here we do for the first time an extensive search for higher overtones of scalar 
perturbations. Some of the lowest lying modes we find are shown in Tables 1, 2 and 3 for 
large, intermediate and small black holes, respectively. 

(i) Large black holes - As proven by Horowitz and Hubeny ^S] in the large black hole 
regime the frequencies must scale as the horizon radius (this can also be proven easily 
and directly from the differential equation (J2J)). We show in Table 1 the results for a 
spherically symmetric mode (7 = 0) for a black hole with r + = 100 which is therefore 
sufficient to infer the behaviour of all large black holes. The fundamental frequency agrees 
with previous results [T^]. Perhaps the most interesting result in this large black hole 
regime is that asymptotically for high overtone number n the frequencies become evenly 
spaced and behave like, for I = 0, 

— = (1.299 - 2.25z)n + 1.856 - 2.673i , (n,r+) ^ oo . (11) 
r+ 

Thus the spacing between frequencies is 

u Sn+1 -u Sn = / 1 299 _ 2 25 -) ( n;r+ )^ 00 . (12) 

Moreover, although the offset 1.856— 2. 673z in (fTTj) is /-dependent (this number is different 
for I = 1 scalar perturbations for example), this asymptotic behaviour for the spacing 
f!12|) holds for any value of Z. In fact our search of the QN frequencies for higher values 
of I reveal that the results are very very similar to those in Table 1. We have gone up to 
I = 4 for scalar perturbations and the results were quite insensitive to I. The asymptotic 
behaviour sets in very quickly as one increases the mode number n. Typically for n = 10 
equation (jll|) already gives a very good approximation. Indeed, for n = 10 we find 
numerically (see Table 1) u s = 1486.23753 — 2516.90740? for a r + = 100 black hole, while 
the asymptotic expression gives u s = 1484.6 — 2517.3i 

(ii) Intermediate black holes - In Table 2 we show some of the lowest lying scalar 
QN frequencies for an intermediate black hole with r + = 1. For a black with this size, 
one finds again that the spacing does not depend on the angular number I for very high 
overtone number n. With an error of about 2% the limiting value for the frequency is, 
for I = 0, 

lu s ~ (1.97- 2.35z> + 2.76- 2.7z , n -> oo . (13) 

For QN frequencies belonging to different Vs the offset in (|T3*j) is different, but as far as 
we can tell numerically, not the asymptotic spacing implied by (|13|). Expression (|T5j) for 
the asymptotic behaviour works well again for n > 10. 

(iii) Small black holes - Our search for the QN frequencies of small black holes, i.e, 
black holes with r + < 1 revealed what was expected on physical grounds, and was 
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Table 1: QN frequencies corresponding to I — scalar perturbations of a large 
Schwarzschild-AdS BH (r+ = 100). It can be seen that for large n the modes become 
evenly spaced. Although not shown here, our numerical data indicates that this happens 
for all values of I and also that the spacing is the same, regardless of the value of /. For 
/ = and for high n the QN frequencies go like f - = (1.299 - 2.25i)n + 1.856 - 2.673i. 
The corresponding spacing between consecutive modes seems to be /-independent. 



n Re[u QN ]: Ioi[ljq N ] : 


n Re[u QN ]: lm[ujQ N ] : 


184.95344 -266.38560 


7 1096.44876 -1841.88813 


1 316.14466 -491.64354 


8 1226.38317 -2066.89596 


2 446.46153 -716.75722 


9 1356.31222 -2291.90222 


3 576.55983 -941.81253 


10 1486.23753 -2516.90740 


4 706.57518 -1166.8440 


50 6682.78814 -11516.9823 


5 836.55136 -1391.8641 


299 39030.810 -67542.308 


6 966.50635 -1616.8779 


300 39160.7272 -67767.3091 



Table 2: QN frequencies corresponding to I = scalar perturbations of an intermediate 
Schwarzschild-AdS BH (r+ = 1). Asymptotically for large n one finds approximately 
u s ~ (1.97 - 2.35i)n + 2.76 - 2.7i. 



n Re[uj QN ]: lra[u QN ] : 


n Re^Qjv]: Im[cj QA r] : 


2.7982 -2.6712 


10 22.44671 -26.20913 


1 4.75849 -5.03757 


11 24.41443 - 28.55989 


2 6.71927 -7.39449 


12 26.38230 -30.91059 


3 8.46153 -9.74852 


13 28.35029 -33.26123 


4 10.6467 -12.1012 


14 30.31839 -35.61183 


5 12.6121 -14.4533 


15 32.28658 -37.96238 


6 14.5782 -16.8049 


16 34.25485 -40.31290 


7 16.5449 -19.1562 


17 36.22318 -42.66340 


8 18.5119 -21.5073 


18 38.19157 -45.01387 


9 20.4792 -23.8583 


19 40.16002 -47.36431 



Table 3: QN frequencies corresponding to / = scalar perturbations of a small 
Schwarzschild-AdS BH (r + = 0.2). Asymptotically for large n one finds approximately 
m s ~ (1.69 - 0.57i)n + 2.29 - OAGi. 



n Re[u QN ]: Im{u QN ] : 


n Re[u QN ]: Im[^ QA r] : 


2.47511 -0.38990 


6 12.45222 -3.89179 


1 4.07086 -0.98966 


7 14.14065 -4.46714 


2 5.72783 -1.57600 


8 15.83026 -5.04186 


3 7.40091 -2.15869 


9 17.52070 -5.61610 


4 9.08118 -2.73809 


10 19.21191 -6.18997 


5 10.7655 -3.31557 


11 20.90359 -6.76355 
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uncovered numerically for the first time in JE] for the fundamental mode: for small black 
holes, the QN frequencies approach the frequencies of pure AdS spacetime [2TJ 123 (see 
also Appendix A). In fact we find 

u s = 2n + I + 3 , r + -> . (14) 

In Table 3 we show some results for a small black hole with r + = 0.2. We stress that the 
values presented in Table 3 for the asymptotic spacing between modes may have an error 
of about 2%. In fact it is extremely difficult to find very high overtones of small black 
holes, and so it is hard to give a precise extimate of the value they asymptote to. 

In summary, we can say that the QN frequencies tend to be evenly spaced asymp- 
totically as n gets very large, no matter if the black hole is large, intermediate or small. 
Moreover the spacing between consecutive modes is, as far as we can tell, independent of 
the angular quantum number I. 



3.2 Electromagnetic quasinormal frequencies 

The fundamental electromagnetic QN frequencies were computed for the first time by 
Cardoso and Lemos [T7j. Recently Berti and Kokkotas JTHJ have redone the calculation 
showing excellent agreement. Here we extend the results to higher overtones. Some of 
the lowest lying electromagnetic frequencies are shown in Tables 4-8. 

(i) Large black holes - As found in ^7] large black holes show a somewhat peculiar 
behaviour: some of the lowest lying modes have pure imaginary frequencies, and these 
are well described by an analytical formula ^7j. A surprising aspect unveiled for the first 
time by the present search is that the number of such modes decreases as the horizon 
radius becomes smaller, as can be seen from Tables 4 and 5. In other words, for very 
large black holes the number of imaginary modes grows. For example, for r + = 1000 
(Table 4) there are eight pure imaginary modes, for r + = 100 there are four such modes 
(see Table 5), and for r + = 10 there are only two. If one wants to go for r+ larger than 
1000, the computation is very time consuming since we use a trial and error method for 
finding new modes. However, not only is this a completely new piece of data, it also 
makes us think that infinitely large black holes may have pure imaginary electromagnetic 
QN frequencies for any overtone number. Perhaps, an infinitly large black hole cannot 
vibrate at all. 

Again, we find that for large black holes and 1 = 1, the frequencies are evenly spaced 
with 

— = (1.299 -2.25z)n- 11.501 + 12i , (n,r + )^oo. (15) 
And a spacing given by 

QJem ° +1 "" cm " =(1.299 - 2.250 , (n,r+)^oo. (16) 
r+ 

For different values of the angular quantum number I, we find the same spacing ()16|) 
between consecutive modes, although the offset in (fTK|) depends on /. So, asymptotically 
for large n and large horizon radius the spacing is the same as for the scalar case! This is 
surprising, specially since the behaviour of the scalar and electromagnetic potentials are 
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Table 4: QNMs corresponding to I = 1 electromagnetic perturbations of a large 
Schwarzschild-AdS BH (r + = 1000). Notice that now there are eight pure imaginary 
modes, still well described by Liu's formula. 



n 


Re[u QN ): 


lm[u QN ] : 


n 


Re[u QN ): 


ImfcuQAr] : 








-1500.004789 


5 





-8985.232 


1 





-2999.982599 


6 





-10596.03 


2 





-4500.093600 


7 





-11644.76 


3 





-5999.513176 


8 


1219.7 


-13566.42 


4 





-7502.69385 


9 


2494.6 


-15847.06 



Table 5: QNMs corresponding to I — 1 electromagnetic perturbations of a large 
Schwarzschild-AdS BH (r + = 100). The first four modes are pure imaginary and are 
well described by Liu's approximation \TJ\. For high n the QN freuencies obey, for 1 = 1, 
i y SSL = (1.299 — 2.25i)n — 11.501 + 12i. The corresponding spacing between consecutive 
modes seems to be /-independent. 



n Re[u QN ): hn[u QN ) : 


n Re[u QN ): lia[u QN ) : 


-150.0479 


10 799.6171 -2171.826 


1 -299.8263 


11 927.812 -2398.208 


2 -450.9458 


12 1056.153 -2624.438 


3 -595.3691 


13 1184.620 -2850.543 


4 22.504 -799.194 


14 1313.192 -3076.546 


5 162.256 -1035.098 


15 1441.856 -3302.464 


6 289.028 -1263.537 


16 1570.601 -3528.310 


7 416.247 -1491.223 


17 1699.416 -3754.094 


8 543.792 -1718.409 


18 1828.295 -3979.824 


9 671.598 -1945.246 


19 1957.229 -4205.508 



radically different. It is even more surprising the fact that this asymptotic behaviour does 
not depend on /, as the electromagnetic potential is strongly /-dependent. Furthermore 
from the first electromagnetic overtones one could surely not anticipate this behaviour. 

(ii) Intermediate black holes - In Table 6 we show some of the lowest lying elec- 
tromagnetic QN frequencies for an intermediate black hole with r + = 1. For a black 
with this size, one finds again that the spacing does not depend on the angular number 
/ for very high overtone number n. With an error of about 2% the limiting value for the 
frequency is, for / = 1, 

u em ~ (1.96 - 2.36i)n + 1.45 - 2.1i , n -»• oo . (17) 

We note that here too the offset in ()17j) does depend on /, but not the asymptotic spacing. 

(iii) Small black holes - For small black holes, see Tables 7 and 8, the spacing seems 
also to be equal as for the scalar case, but since it is very difficult to go very high in mode 
number n in this regime, the error associated in estimating the asymptotic behaviour is 
higher, and one cannot be completely sure. Again, the electromagnetic QN frequencies of 
very small black holes asymptote to the pure AdS electromagnetic modes (see Appendix 
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Table 6: QNMs corresponding to I — 1 electromagnetic perturbations of an intermediate 
Schwarzschild-AdS BH (r+ = 1). Asymptotically for large n the modes become evenly 
spaced in mode number and behave as cu em ~ (1.96 — 2.36i)n + 1.45 — 2.H. 



n ne[ujQN\. im[ujQiy\ . 




2.163023 -1.699093 


10 21.067466 -25.61714 


1 3.843819 -4.151936 


11 23.015470 -27.98278 


2 5.673473 -6.576456 


12 24.965381 -30.34713 


3 7.553724 -8.980538 


13 26.916889 -32.71037 


4 9.458385 -11.37238 


14 28.869756 -35.07265 


5 11.37722 -13.75633 


15 30.823790 -37.43413 


6 13.30526 -16.13482 


16 32.778838 -39.79488 


7 15.23974 -18.50933 


17 34.734776 -42.15499 


8 17.17894 -20.88081 


18 36.691500 -44.51455 


9 19.12177 -23.24993 


19 38.648922 -46.87360 



Table 7: QNMs corresponding to / = 1 electromagnetic perturbations of a small 
Schwarzschild-AdS BH (r + = 0.2). Asymptotically for large n one finds approximately 
u em ~ (1.68 - 0.59i)n + 1.87 - O.OAi. 



n Re[uj QN ]: lra[u QN ] : 


n He [ujqn] '■ lra[uJ QN ] : 


2.63842 -0.05795 


6 12.00066 -3.53148 


1 3.99070 -0.47770 


7 13.66436 -4.12974 


2 5.49193 -1.08951 


8 15.33370 -4.72479 


3 7.07835 -1.70859 


9 17.00715 -5.31725 


4 8.70165 -2.32191 


10 18.68370 -5.90758 


5 10.3450 -2.92920 


11 20.36268 -6.49615 



Table 8: The fundamental (n = 0) QNMs corresponding to I — 1 electromagnetic 
perturbations of a small Schwarzschild-AdS BH for several values of r+. 



r+ Re[uQ N ]: Im^Qjv] : 


r+ Re[u QN ]: lm[u QN } : 


1/2 2.25913 -0.65731 


1/10 2.85188 -0.00064 


1/3 2.40171 -0.29814 


1/12 2.88058 -0.00030 


1/4 2.53362 -0.13364 


1/16 2.91363 -0.00016 


1/5 2.63842 -0.05795 


1/18 2.92406 -0.00009 


1/8 2.80442 -0.00565 


1/20 2.93200 -0.00002 
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A, where we sketch their computation). Indeed we find that 

u emAdS = 2n + l + 2 ,r + -> 0. (18) 

This can be clearly seen from Table 8, where we show the fundamental mode for small 
black holes of decreasing radius. As the horizon radius gets smaller and smaller, the 
fundamental frequency approaches the value of 3 + Oi, which is indeed the correct pure 
AdS mode for I = 1, n = 0, electromagnetic perturbations. It was conjectured by 
Horowitz and Hubeny 15 a that for very small black holes in AdS space, the imaginary 
part of the QN frequency for spherically symmetric perturbations should scale with the 
horizon area, i.e., with ri. Their argument was based on a previous result for the 
absorption cross section for the I = component. This conjecture was later verified 
numerically to be correct by Konoplya ^Hj for the / = case. ^From Table 8 it is 
however apparent that this scaling is no longer valid for I = 1 perturbation, and indeed 
we find it is not valid for I ^ perturbations, be it scalar, electromagnetic or gravitational 
perturbations. The reason why the imaginary part no longer scales with the horizon area 
for I ^ perturbations is due to the fact that the partial absorption cross section only 
scales with the horizon area for I = perturbations. For other Z's the behaviour is more 
complex, and it could be that there is no simple scaling, or even that the behaviour is 
oscillatory with the mass M of the black hole. We refer the reader to j21] for details on 
the absorption cross section of black holes. 



3.3 Gravitational quasinormal frequencies 

The fundamental gravitational QN frequencies were computed for the first time by Car- 
doso and Lemos [TZj. We remind that there are two sets of gravitational wave equations, 
the odd and even ones. Although it was found [Tj\ that there is a family of the odd 
modes which is very slowly damped and pure imaginary, it was possible to prove that for 
high frequencies both odd and even perturbations must yield the same QN frequencies. 
We present the results for higher overtones of odd perturbations in Tables 9-12, and even 
perturbations in Tables 13-16. 

3.3.1 Odd perturbations 

(i) Large black holes - As discussed for the first time in [T2| these exhibit a pure 
imaginary fundamental mode (see Table 9.). For large black holes, this mode is slowly 
damped and scales as the inverse of the horizon radius. Our analysis for higher Z's 
indicates that in the large black hole regime an excellent fit to this fundamental pure 
imaginary mode is 

(l-l)(l + 2) , , 

^odd^o = 1 > r + -> 00 • ( 19 ) 

This generalizes a previous result by Berti and Kokkotas fH] for the I = 2 case. The 
simplicity of this formula (which is just a fit to our numerical data), leads us to believe 
it is possible to find an analytical explanation for it, but such expalnation is still lacking. 
In the large black hole regime, asymptotically for high overtones one finds, for / = 2 for 
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example, 

— = (1.299 -2.25i)n + 0.58- 0.421 , (n,r+)^oo. (20) 
r+ 

This leads to the spacing 

UJoddn+1 ~ Uoddn = (1.299 - 2.25i) , (n,r + )^oo, (21) 
r + 

which, as our results indicate is again /-independent. Again, the offset in (J2U|l depends 
on Z. 

(ii) Intermediate black holes - Results for the odd QN frequencies of an intermediate 
(r + = 1) black hole are shown in Table 10. With an error of about 5% the limiting value 
for the frequency is, for 1 = 2, 

u dd ~ (1-97 - 2.35i)n + 0.93 - 0.32z , n -> oo . (22) 



We note that here too the offset in ()22|) does depend on I, but not the asymptotic spacing, 
with a numerical error of about 5%. 

(iii) Small black holes - The behavior for small black holes is shown in Tables 11 and 
12. As the black hole gets smaller, the pure imaginary mode gets more damped: the 
imaginary part increases, as can be seen from Table 12, where we show the two lowest 
QN frequencies for small black holes with decreasing radius. As mentioned by Berti and 
Kokkotas JB] the ordering of the modes here should be different. However, since one can 
clearly distinguish this pure imaginary mode as belonging to a special family, we shall 
continue to label it with n = 0. We have not been able to follow this mode for black holes 
with r + < 0.5, and so Table 12 does not show any pure imaginary modes for horizon 
radius smaller than 0.5. We note that, as for the scalar and electromagnetic cases, here 
too the modes are evenly spaced, with a spacing which seems to be independent of I no 
matter if the black hole is large or small. For very small black holes, the frequencies 
reduce to their pure AdS values, computed in Appendix A, to wit 

tu odd = 2n + l + 2, (23) 

One can see this more clearly from Table 12, where in fact for very small black holes the 
frequency rapidly approaches (J23|) . Again, for small black holes, the imaginary part does 
not scale with the horizon area, by the reasons explained before, in section 3.2. 

In conclusion, the higher overtones of odd perturbations follow a pattern very similar 
to the scalar case. We note that the asymptotic behaviour sets in very quickly, much 
like what happened for scalar and electromagnetic perturbations. Typically the formulas 
yielding the asymptotic behaviour work quite well for n > 10. We are now able to 
prove that for sufficiently high frequencies the scalar and gravitational perturbations are 
isospectral, a mystery that remained in [T7], This is done in section 4.1 below. 

3.3.2 Even perturbations 

Let us now briefly discuss the even modes. As found previously [T7| these modes behave 
very similar to the scalar ones. Yet, the even gravitational modes are stipulated by a 
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Table 9: QNMs corresponding to I = 2 odd gravitational perturbations of a large 
Schwarzschild-AdS BH (r + = 100). The fundamental QN frequency is pure imaginary 
and seems to be well described by the formula uj n=0 = ~ (7 ~ 1 - >( ^ +2 - > i valid only in the large 
black hole regime. In the large n limit one finds ^ = (1.299 - 2.25i)n + 0.58 - 0.42i 
The corresponding spacing between consecutive modes seems to be /-independent. 



n Re[u QN ]: lm[ujQ N } : 


n Re[u QN }: lm[u QN ] : 


- 0.013255 


6 836.55392 -1391.86345 


1 184.95898 - 266.38403 


7 966.50872 -1616.87735 


2 316.14887 - 491.64242 


8 1096.45098 -1841.88755 


3 446.46505 - 716.75629 


9 1226.38527 -2066.89540 


4 576.56293 - 941.81172 


10 1356.31422 -2291.90170 


5 706.57797 - 1166.8433 


50 6552.87704 -11291.9807 



Table 10: QNMs corresponding to I — 2 odd gravitational perturbations of an intermedi- 
ate Schwarzschild-AdS BH (r+ = 1). Asymptotically for large n one finds approximately 
^odd ~ (1-97 - 2.35i)n + 0.93 - 0.32i. 



n Re[u QN ]: lm[u QN ] : 


n Re[u QN }: Iui[ujq N ] : 


-2 


10 20.604949 -23.803860 


1 3.033114 -2.404234 


11 22.567854 -26.157246 


2 4.960729 -4.898194 


12 24.531429 -28.510214 


3 6.905358 -7.289727 


13 26.495564 -30.862849 


4 8.854700 -9.660424 


14 28.460169 -33.215214 


5 10.80784 -12.02344 


15 30.425175 -35.567355 


6 12.76384 -14.38266 


16 32.390524 -37.919308 


7 14.72199 -16.73969 


17 34.356173 -40.271103 


8 16.68179 -19.09530 


18 36.322082 -42.622761 


9 18.64286 -21.44994 


19 38.288221 -44.974301 



Table 11: QNMs corresponding to / = 2 odd gravitational perturbations of a small 
Schwarzschild-AdS BH (r + = 0.2). Asymptotically for large n one finds approximately 
o^odd - (1-69 - 0.59i)n + 2.49 + 0.06i. 



n Re[u QN ]: lm[u QN ] : 


n Re[iVQ N ]: lm[ujQ N } : 


2.404 -3.033 


6 12.67161 -3.43609 


1 4.91594 -0.30408 


7 14.33020 -4.05366 


2 6.30329 -0.89773 


8 15.99881 -4.66448 


3 7.82330 -1.53726 


9 17.67433 -5.26955 


4 9.40720 -2.17744 


10 19.35465 -5.86978 


5 11.0279 -2.81083 


11 21.03839 -6.46596 
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Table 12: The fundamental (n = 0) QNMs corresponding to I = 2 odd gravitational 
perturbations of a small Schwarzschild-AdS BH for several values of r + . 



r+ 




Re[u QN \: 


Im[u;QAr] : 


r+ 


Re[u QN ]: 


ha[uQ N ] : 


0.8 (n 


= 0) 





-3.045373 


0.5 (n = 1) 


3.03759 


-0.71818 


0.8 (n 


= 1) 


2.89739 


-1.69556 


0.4 


3.16209 


-0.43092 


0.7 (n 


= 0) 





-3.83538 


0.3 


3.35487 


-0.17320 


0.7 (n 


= 1) 


2.90665 


-1.34656 


0.2 


3.62697 


-0.01792 


0.6 (n 


= 0) 





-4.901973 


0.1 


3.84839 


-0.00005 


0.6 (n 


= 1) 


2.95550 


-1.02196 


1/15 


3.90328 


-0.00001 


0.5 (n 


= 0) 





-6.40000 


1/20 


3.92882 


-0.000002 



more complicated potential, and we have to truncate the series in power of x — x + at 
larger N, which makes the whole procedure be more time consuming. That is why when 
considering small black holes we were restricted only by first seven modes in that case. 
It is, however, sufficient to see that even gravitational QNMs, similar to other kind of 
perturbations, tend to arrange into equidistant spectrum under the increasing of n. We 
show in Tables 13-16 the numerical results for the QN frequencies of even gravitational 
perturbations. 

(i) Large black holes - Results for the QN frequencies of large black holes are shown 
in Table 13. In this regime one finds for I = 2 even perturbations 

= (1.299 -2.25i)n + 1.88 -2.661 , (n,r+) ^ oo , (24) 

r + 

leading to the spacing 

u evean+1 - ow n = (L299 _ 2>25j) ? ( n ,r + )^oo, (25) 
r + 

which once more turns out to be /-independent! All the results concerning the spacing 
of frequencies for large black holes have a very good precision, since in this regime it is 
possible to go very far out in overtone number (typically n = 300 is enough to achieve a 
0.1% accuracy for the spacing). 

(ii) Intermediate black holes - In Table 14 we show some of the lowest lying even 
gravitational QN frequencies for an intermediate black hole with r + = 1. For a black 
with this size, one finds again that the spacing does not seem to depend on the angular 
number I for very high overtone number n. With an error of about 5% the limiting value 
for the frequency is, for 1 = 2, 

u even ~ (1.96 - 2.35i)n + 2.01- 1.5i , n -> oo . (26) 

We note that here too the offset in (fTTjl does depend on Z, but not the asymptotic spacing. 

(iii) Small black holes - The behavior for small black holes is shown in Tables 15 and 
16. Our search for the QN frequencies of small black holes, i.e, black holes with r + < 1 
revealed again what was expected on physical grounds: for small black holes, the QN 
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Table 13: QNMs corresponding to I — 2 even gravitational perturbations of a large 
Schwarzschild-AdS BH (r+ = 100). For large n, one finds ^ = (1.299 -2.25i)ra + 0.58- 
0.42i. The corresponding spacing between consecutive modes seems to be /-independent. 



n Re[u;Qjv]: Im[cjQjv] : 


n Re[u QN ]: ha[uQ N ] : 


184.97400 -266.351393 


6 966.609780 -1616.695872 


1 316.17838 -491.584999 


7 1096.56635 -1841.681256 


2 446.50884 -716.674054 


8 1226.51495 -2066.664293 


3 576.62103 -941.70468 


9 1356.45821 -2291.645761 


4 706.65039 -1166.71147 


10 1486.39776 -2516.626168 


5 836.64066 -1391.70679 


50 6683.51993 -11515.70869 



Table 14: QNMs corresponding to I — 2 even gravitational perturbations of an intermedi- 
ate Schwarzschild-AdS BH (r+ = 1). Asymptotically for large n one finds approximately 
eleven ~ (1-96 - 2.35i)n + 2.01 - L5i 



n Re[u QN ]: lra[uJ QN } : 


n Re[u QN ): lm[u QN ] : 


3.017795 -1.583879 


10 21.68949 -24.98271 


1 4.559333 -3.810220 


11 23.64402 -27.33549 


2 6.318337 -6.146587 


12 25.60052 -29.68799 


3 8.168524 -8.500194 


13 27.55860 -32.04026 


4 10.061220 -10.85631 


14 29.51796 -34.39234 


5 11.976813 -13.21224 


15 31.47838 -36.74424 


6 13.906140 -15.56749 


16 33.43969 -39.09600 


7 15.844371 -17.92208 


17 35.40174 -41.44762 


8 17.788721 -20.27609 


18 37.36444 -43.79914 


9 19.737469 -22.62960 


19 39.32769 -46.15057 



Table 15: QNMs corresponding to I = 2 even gravitational perturbations of a small 
Schwarzschild-AdS BH (r + = 0.2). Asymptotically for large n one finds approximately 
eleven ~ (1-61 - O.Qi)n + 2.7 + 0.37i. 



n 


Re[u QN ]: 


lm[u QN ] : 


n 


Re[u QN ]: 


Im[^ QA r] : 





3.56571 


-0.01432 


3 


7.65872 


-1.42994 


1 


4.83170 


-0.26470 


4 


9.20424 


-2.04345 


2 


6.17832 


-0.82063 


5 


10.78800 


-2.65360 



Table 16: The fundamental (n = 0) QNMs corresponding to I — 2 even gravitational 
perturbations of a small Schwarzschild-AdS BH for several values of r+. 



r + Bb[u qn ]: lra[u QN ] : 


r + Re[u QN ]: lra[ujQ N ] : 


0.8 2.91541 -1.18894 


0.3 3.29299 -0.14103 


0.7 2.90591 -0.98953 


0.2 3.56571 -0.01432 


0.6 2.92854 -0.78438 


0.1 3.80611 -0.00005 


0.5 2.98985 -0.57089 


1/15 3.8735 -0.00001 


0.4 3.10317 -0.35043 


1/20 3.90852 -0.000002 
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frequencies approach the frequencies of pure AdS spacetime (see Appendix A). In fact 
we find 

^cven AdS = 2n + I + 2 , r+ -> . (27) 

In Table 15 we show the lowest lying QN frequencies for a small black hole (r+ = 0.2). We 
stress that the values presented in Table 15 (as a matter of fact, all the Tables containing 
data for small black holes) for the asymptotic spacing between modes may have an error 
of about 2%. In fact it is extremely difficult to find very high overtones of small black 
holes, and so it is hard to give a precise extimate of the value they asymptote to. In 
Table 16 we show some the fundamental even QN frequencies for small black holes of 
decreasing radius, and one can clearly see how the fundamental frequency approaches the 
pure AdS value given in Appendix A. 



4 Discussion of the results 



4.1 Why are the scalar and gravitational perturbations isospec- 
tral in the large black hole regime? 

In a previous paper (section IIIC in [E|), we have showed why the odd and even gravita- 
tional perturbations yield the same QN frequencies for large frequencies. The whole ap- 
proach was based on the fact that the odd and even gravitational potentials are superpart- 
ner potentials i.e., they are related to one another via V Q( id = W 2 +j^- +0 , V even = 
W 2 -f^+f3 , where (3 = The function W is W = ffi+^+ ^g+ffi 3 - 

I (H + + ^a~)- ^ or more details we refer the reader to [T7]. We shall now see that 
a similar method can be applied to show that in the large black hole regime, scalar and 
gravitational perturbations are isospectral for large QN frequencies. To begin with, we 
note that the potentials VI and VI defined by 



and 



. a 6M\ 



with / = jjff+f — — , and a any constant, are superpartner potentials. The superpotential 
W is in this case is given by 

Thus, the two superpartner potentials VI and V2 can be expressed in terms of W as 

V1 = W 2 +—, V2 = W 2 — — . 31 

Why are these two potentials of any interest? Because in the large r + limit, which we 
shall take to be r + ^> a, we have / ~ r 2 — — . Notice now that in this large r + limit 
the scalar potential (J3J) is V & ~ f(2 + 2 ^r), with / ~ r 2 — — , since in this limit and 
with r + 3> I, one has l( ~ l+ 2 1 ^ 2. Thus, VI reduces to the scalar potential and V2 to 
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the gravitational odd potential, provided we take a = 1(1 + 1). It then follows from the 
analysis in [T7] (section IIIC) that for large black holes these two potentials should yield 
the same frequencies. 

4.2 Future Directions 

The preceding sections have shown that the QNMs of Schwarzschild-AdS black holes 
have a universal behavior in the asymptotic regime of high overtones. This was verified 
explicitly and with great accuracy for the large black hole regime, where we showed 
numerically that the spacing does not depend on the perturbation in question and is 
equal to 

Un+1 ~ Un ' = (1.299 - 2.25i) , (n,r+) ^ oo . (32) 

We conjecture that the asymptotic behavior is the same for all kinds of perturbations 
irrespectively of the black hole size, i.e., a fixed horizon radius r + Schwarzschild-AdS 
black hole will have an asymptotic spacing between consecutive QN frequencies which is 
the same for scalar, electromagnetic and gravitational perturbations. The difficulty in 
extracting very high overtones for small black holes however, prevents us from having an 
irrefutable numerical proof of this. It would be extremely valuable to have some kind 
of analytical scheme for extracting the asymptotic behavior, much as has been done for 
the asymptotically flat space by Motl and Neitzke |3]. However it looks quite difficult to 
make any analytical approximation in asymptotically AdS spaces, although there have 
been some atempts at this recently (see for example Musiri and Siopsis [Sj). We also 
note that the spacing (}3*2*|) was already found to be true by Berti and Kokkotas [THJ for 
the scalar and gravitational cases for the lowest radiatable multipole, i.e, I = and I = 2 
scalar and gravitational perturbations respectively. We have concluded that, surprisingly, 
the spacing also works for electromagnetic case and for any value of I. It was 
observed that, despite having such different potentials the scalar, the electromagnetic and 
gravitational QN frequencies have the same asymptotic behavior. Can one formulate some 
very general conditions the potentials should obey in order to have the same asymptotic 
solutions? This is still an open question. 

There has been recently an exciting development trying to relate the asymptotic 
QN frequencies with the Barbero-Immirzi parameter [HE]. In fact it was observed, in 
the Schwarzschild case, that asymptotically for high overtones, the real part of the QN 
frequencies was a constant, /-independent, and using some (not very clear yet) correspon- 
dence between classical and quantum states, was just the right constant to make Loop 
Quantum Gravity give the correct result for the black hole entropy. Of course it is only 
natural to ask whether such kind of numerical coincidence holds for other spacetimes. 
We have seen that apparently we are facing, in AdS space, a universal behavior, i.e., the 
asymptotic QN frequencies do not depend on the kind of perturbations, and also don't 
depend on I. However, and in contrast with asymptotically flat space, the real part of the 
asymptotic QN frequency is not a constant, but rather increases linearly with the mode 
number n. This is no reason to throw off the initial motivation of seeking some kind of 
relation between Loop Quantum Gravity and QNMs, after all, there are no predictions 
for AdS space. 
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Finally we point out that the asymptotic behavior studied here for the Schwarzschild- 
AdS black hole will hold also for other black holes in asymptotically AdS. One example 
of these is the black hole with non-trivial topology jTHj . The general line element for this 
spacetime is [TU] : 

ds 2 = f{r) dt 2 - fir^dr 2 - r 2 (dO 2 + dcj) 2 ) (33) 

where 

„ , r 2 AMR . ^ 

f u = *—r ' (34) 

where M is the ADM mass of the black hole, and R is the AdS radius. There is a 
horizon at r + = (4M) 1 / 3 i?. The range of the coordinates 9 and (ft dictates the topology 
of the black hole spacetime. For a black hole with toroidal topology, a toroidal black 
hole, the coordinate 9 ranges from to 2ir, and <fi ranges from to 2tt as well. For the 
cylindrical black hole, or black string, the coordinate 9 has range —oo<R6<oo, and 
< < 2tt. For the planar black hole, or black membrane, the coordinate is further 
decompactified — oo < R<p < oo jTHj. The fundamental QN frequencies for these black 
holes were computed in |2H|, where it was verified that they follow the same pattern as 
for Schwarzschild-AdS black holes. Indeed one easily sees that in the large black hole 
regime they both should yield the same results as the potentials are equal in this regime 
(compare the potentials in [20] with the ones in the present work). In particular the 
asymptotic behavior will be the same. 



5 Conclusion 

We have done an extensive search for higher overtones n of the QNMs of Schwarzschild- 
AdS BH corresponding to scalar, electromagnetic, and gravitational perturbations. We 
have shown that: (i) No matter what the size of the black hole is, the QN frequencies 
are evenly spaced, both in the real and in the imaginary component, for high overtone 
number n; (ii) The spacing between consecutive modes is independent of the perturba- 
tion. This means that scalar, electromagnetic and gravitational perturbations all have, 
asymptotically, the same spacing between modes. This is one of the major findings in this 
work, together with the fact that this spacing seems to be also independent of the angu- 
lar quantum number I; (iii) We were able to prove that the scalar and gravitational QN 
frequencies must asymptotically be the same; (iv) The electromagnetic QN frequencies 
of large black holes have a number of first overtones with pure imaginary parts, and the 
higher the black hole radius r + , the higher the number of these first pure damped, non- 
oscillating modes; (v) Finally, we have computed analytically the electromagnetic and 
gravitational pure AdS modes, and we have shown numerically that the QN frequencies 
of very small black holes asymptote to these pure AdS modes; 
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A Pure AdS Normal Modes for Electromagnetic and 
Gravitational Perturbations 

In this appendix we shall briefly outline how to compute the pure modes of AdS space 
(no black hole, M — 0) for electromagnetic and gravitational perturbations. The scalar 
case was dealt with by Burgess and Lutken [21] (see also [22])- In pure AdS space the 
electromagnetic and gravitational potentials (both odd and even) are 

as can be seen by substituting M = in (JSJ)- (JI|). Also in this case the relation r(r*) 
takes the simple form 

T 

r = Rtan — , (36) 
R 

and therefore the potential (}35|) takes a simple form in the r* coordinate, namely 

V = -?£±V (37) 

To proceed, we note that the change of variable x = sin (j^j leads the wave equation to 
a hypergeometric equation, 

+ ^ + -^) = 0, (38) 

ox z a ox a z 

with 

a = A(ujRfx(l -x)- 4:1(1 + 1)(1 - x) , (39) 

a = Ax(l-x), (40) 
f = 2(l-2x). (41) 

To put this in a more standard form, one changes wavefunction by defining 

f(x) = Va^I# Z(x), (42) 
and one gets the following standard hypergeometric differential equation for Z: 

d 2 Z(x) dZix) . , 
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with a defined in (J40|) and 



r = 6 -4/(x - 1) - I2x, (44) 
A = -4 - 4/ - I 2 + u 2 . (45) 

By requiring well behaved fields everywhere a simple analysis j2H] then shows that the 
following constraint needs to be satisfied, 

uj R = 2n + I + 2 . (46) 

These are the pure AdS frequencies for electromagnetic and gravitational perturbations, 
corresponding to pure AdS normal modes of the corresponding fields. One can compare 
the frequencies in (I46|) with the scalar frequencies corresponding to pure AdS modes 
HUE2I, oj s R = 2n + Z + 3. 
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